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Abstract 



We take a step toward a nonperturbative gravitational path integral for black- 
hole geometries by deriving an expression for the expansion rate of null geodesic 
^ i congruences in the approach of causal dynamical triangulations. We propose to 

I use the integrated expansion rate in building a quantum horizon finder in the sum 

over spacetime geometries. It takes the form of a counting formula for various 
types of discrete building blocks which differ in how they focus and defocus light 
rays. In the course of the derivation, we introduce the concept of a Lorentzian 
dynamical triangulation of product type, whose applicability goes beyond that of 
describing black-hole configurations. 
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1. Introduction 

Very encouraging progress has been made recently in constructing spacetime dy- 
namically from a nonperturbative gravitational path integral, by studying the 
continuum limit of causal dynamical triangulations ^1 12 El El- The quantum 
geometries generated in this way exhibit semiclassical properties at sufficiently 
large scales: they are four-dimensional E] and the large-scale dynamics of their 
spatial volume is described by an effective cosmological minisuperspace action 
[?]. Their short-distance behaviour is highly nonclassical, including a smooth 
dynamical reduction of the spectral dimension from four to two Pj and evidence 
of fractality P . 

A question frequently asked of this and other approaches to quantum gravity 
is what they have to say about the quantum dynamics of black holes and, more 
specifically, whether the black-hole entropy "comes out right". What is usually 
meant by this is whether the theory can reproduce the Bekenstein-Hawking for- 
mula S = A/A which relates the entropy 5 of a black hole to its area A in Planck 
units, and preferably at the same time provide a microscopic explanation of the 
presence of black-hole entropy in terms of fundamental excitations of geometry 
(see jniUnHH] for recent reviews). Although often portrayed as a touchstone for 
quantum gravity, it should be kept in mind that the entropy formula and other 
thermodynamic relations satisfied by black holes are semiclassical. Whether or 
not they have a fundamental role to play in a genuinely nonperturbative formu- 
lation of the theory remains to be seen. 

This naturally raises the question whether the approach of causal dynami- 
cal triangulations, which we believe is a strong candidate for a nonperturbative 
theory of four-dimensional quantum gravity, can provide new insights into the 
quantum properties of black holes. Apart from having already reproduced cer- 
tain classical aspects of general relativity from first principles, the fact that the 
causal, Lorentzian structure of spacetime geometry plays a central role and that 
a well-defined Wick rotation is available at the regularized level seem to make 
this formulation particularly suited for addressing issues to do with black holes. 
However, putting this into practice turns out to be a challenging proposition. Not 
only are there technical obstacles to be overcome, but one also has to decide what 
precise quantity should be calculated if one were given a well-defined method for 
performing nonperturbative sums over geometries, such as causal dynamical tri- 
angulations, a question that has hardly begun to be addressed. 

The main choices to be made when setting up the path integral are that of an 
ensemble of geometries to be summed over, and of boundary conditions for the 
geometries. Since numerical evaluations of the path integral must necessarily take 
place in a finite spacetime volume, boundary conditions will have to be provided 
not only on some initial and final spatial slice, but also on a (timelike) boundary 



2 



at large radius, and potentially also on a hypersurface at small radius to avoid 
any central singularity. ^ We will for the moment set aside the question of how the 
boundary conditions should be chosen, and how many black holes one can expect 
to generate dynamically as a function of the boundary data. Instead we will 
focus on an issue that will be relevant regardless of the geometric ensemble and 
boundary conditions chosen, namely, which observable in the quantum theory can 
give us information about the presence or otherwise of a black-hole configuration. 
In a fully nonperturbative formulation, this is a difficult task because of the 
absence of a classical background structure. In the nonperturbative path integral, 
all possible geometries are superposed, not just those which represent fluctuations 
around a given classical background geometry. Semiclassical geometry emerges 
only in the continuum limit and at sufficiently large scales [3 E] . The problem is 
then how and where in the dynamically generated quantum geometry one should 
look for evidence of a black hole, say, an apparent horizon. 

To simplify matters slightly, we will consider geometries with at most a single 
black hole, and which moreover have (an approximate) spherical symmetry. To 
translate these conditions to the setting of causal dynamical triangulations, a use- 
ful concept is that of triangulations of product type. These are roughly speaking 
simplicial analogues of fibred spaces, with a triangulated base space and triangu- 
lated fibres. The inspiration for such structures comes in part from Lorentzian 
semi-random lattices ^21; which in turn were motivated by causal triangulated 
models of quantum gravity in lower dimensions 

In the context of a nonperturbative description of four-dimensional black 
holes, we propose to use triangulations of product type whose base space is a 
simplicial version of the r-t-plane, with fibres representing the two angular direc- 
tions. For simplicity, it seems a good strategy not to include completely general 
geometries (triangulations) of this type in the path integral initially, but only 
a subset which satisfies certain homogeneity requirements along the angular di- 
rections, implementing an approximate spherical symmetry at a coarse-grained 
scale. The fact that an exact continuous rotation symmetry cannot even in prin- 
ciple be realized by the simplicial structures we are working with is not at all 
a drawback, but simply implies that the path integral will necessarily include 
fiuctuations which violate spherical symmetry, especially at short distances. This 
is desirable from the point of view of the quantum theory, because it is likely 

""^A detailed discussion of the inclusion of boundary terms in causal dynamical triangulations 
can be found in [T^ . 

^The presence of a Mobius inversion formula for such semi-random lattices, relating their 
partition function to that of a statistical model of geometric objects of one dimension less (see 
also 53), may help to solve these models analytically. In the context of three-dimensional 
quantum gravity in terms of causal dynamical triangulations, this possibility is currently being 
explored 15 . A related quantum-cosmological model which introduces more "order" in three- 
dimensional random triangulations is described in |16|. 
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to make the model closer to the full theory than gravitational models where an 
exact symmetry reduction is performed before the quantization (see JT7] for an 
analogous reasoning in canonical quantum gravity). 

In the present work, we will not attempt to define and evaluate a nonpertur- 
bative path integral for black holes, but pursue a more modest goal, namely, to 
formulate a geometric observable for dynamical triangulations to help to deter- 
mine the "presence" (in the sense of expectation values) of a black hole in the 
quantum theory. The observable is a simplicial version of the integrated expansion 
rate for light rays, whose vanishing in the classical continuum theory is an indica- 
tor of the presence of an apparent horizon. We derive an analogous quantity for 
causal dynamical triangulations which has a particularly simple form, and which 
we hope can be used to construct an efficient "horizon finder" in Monte Carlo 
simulations of the quantum theory. It depends only on the numbers of simplicial 
building blocks of various types and orientations occurring in a given fibre of the 
product triangulation, but not on how these building blocks are glued together 
locally. In this sense, the presence of an apparent horizon can be established 
simply by counting.^ 

The rest of the paper is structured as follows. In Sec. 2 we introduce the 
general concept of a Lorentzian triangulation of product type, before specializing 
to the case of a two-dimensional base space. For product triangulations in 2 + 1 
and 2 + 2 dimensions, we classify all simplices occurring in the fibres according 
to their orientation and derive some topological relations among them. Sec. 
3 is devoted to a discussion of extrinsic curvatures and expansion rates. We 
start by recalling the classical notion of a trapped surface and compute the light 
expansion rates on a surface of codimension 2 from the extrinsic curvatures and 
other geometric data. We then repeat the calculation for piecewise fiat manifolds, 
by employing a careful limiting process at the "kinks" of the triangulation.^ This 
sets the stage for computing the integrated expansion rates for causal dynamical 
triangulations. We derive the "counting formulas", our main result, in both 2+1 
and 2 + 2 dimensions, and also give a qualitative interpretation in terms of the 
focussing and defocussing of light rays by particular types of simplices. In Sec. 4, 
we construct an explicit example of a triangulated black hole in the formalism of 
Lorentzian product triangulations, and Sec. 5 contains a summary and outlook. 
The appendix deals with the definition of affine coordinates, which are used in 
some of our derivations. - Part of the work presented here is contained in the 
diploma thesis of one of the authors [19j. 

•^To avoid any misunderstandings, let us emphasize that this counting has nothing to do 
with a counting of microstates to obtain an entropy for a black hole. 

''A related treatment in the context of classical Regge calculus has appeared previously in 

m 
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2. Lorentzian triangulations of product type 

A triangulated manifold T of product type is a particular case of a simplicial 
manifold. Topologically speaking, it is a cartesian product T = BxFofab- 
dimensional base space B with an /-dimensional fibre space F. A simplicial 
realization of this structure consists of a simplicial base manifold B, where to 
each fc-dimensional simplex, < k < b, (fc-simplex for short) o" of i? we associate 
a so-called a-tower, which is a particular triangulation of cr x F, whose top- 
dimensional simplices are of dimension k+f. We will only consider triangulations 
T with a finite number of (i-simplices. For k = b, the dimension of the simplices 
in the a-towers is maximal, and coincides with the dimension d = b+f of T. The 
triangulations of these towers as well as the number of (i-simplices they contain 
will in general depend on the 6-dimensional base simplex a. In order that the 
cr-towers fit together in a simplicial manifold, neighbouring towers have to satisfy 
matching conditions along their common (d — l)-dimensional boundaries, which 
themselves are a-towers over the (6— l)-dimensional simplices of the base manifold 
B. The triangulations of the ci-dimensional towers induce triangulations on all 
lower- dimensional a-towers, that is, on the vertex towers, the edge towers, etc. 
By definition, each vertex of a product triangulation lies in precisely one vertex 
tower. 

The focus of our interest will be on Lorentzian triangulations of product type, 
and our main applications will have a two-dimensional base triangulation. Note 
that the simplicial manifolds one considers in the approach of causal dynamical 
triangulations can be thought of as a special class of triangulations of product 
type, with base space B the one-dimensional proper-time direction, and fibre 
F a spatial hypermanifold of constant time. The simplicial realization of B is 
simply a one-dimensional chain of time-like edges, all of equal length. For d- 
dimensional simplicial spacetimes, the spatial slices of constant integer time t are 
the ((i — l)-dimensional vertex towers over the (zero-dimensional) vertices of B, 
and the minimal spacetime "sandwiches" between times t and t + 1 are the edge 
towers over B. Each vertex of the spacetime triangulation lies in the vertex tower 
of exactly one vertex in B and thus inherits an integer time coordinate. 

The simplest non-trivial case of a causal triangulation is in dimension two, an 
example of which is depicted in Fig. ^ The vertex towers are one-dimensional 
chains of spatial edges or links, and the edge towers are linear sequences of two- 
dimensional triangles pointing up or down. In a useful notation that generalizes 
to more complicated examples, we denote an up-triangle by [2, 1] ("two vertices in 
the tower above the base vertex at time t, one vertex in the tower above the base 
vertex at ), and a down-triangle by [1, 2] (the converse). Each strip is uniquely 
characterized by a sequence of such number pairs. Analogously, the ci-simplices 
in a minimal spacetime sandwich of a ^-dimensional causal triangulation T can 
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Figure 1: Example of a two-dimensional Lorentzian triangulation T whose base space 
B is the one-dimensional triangulation consisting of three timelike edges. Note that in 
our graphical representation lengths are not represented isometrically, in fact, this is 
impossible because of the intrinsic curvature carried by T. There is only one type of 
triangular building block which appears with either up- or down-orientation. Conse- 
quently, all spacelike edges (horizontal) and all timelike edges (interpolating between 
adjacent slices of constant integer time) have identical length, although the latter fact 
is not rendered faithfully in the figure. 

be characterized by a pair [ii,i2], h £ {^,---,d}, where ii counts the number 
of vertices at time t and 22 those at time t + 1 4J. Since a (i-simplex has d + 1 
vertices, there are d different ways how its vertices can be distributed over the 
two constant-time sUces, leading to configurations of type [d, 1], [d—1, 2], [l,d\. 

The product triangulations we will consider are Lorentzian, that is, they are 
assembled from flat d-dimensional Minkowskian simplicial building blocks and 
have a product structure with respect to the time direction. We show in the 
appendix that in such manifolds the notion of a constant-time slice can be ex- 
tended naturally to non-integer t. The resulting (rf— l)-dimensional hypersurfaces 
are again piecewise flat, but the individual building blocks are not necessar- 
ily simplices (for example, the intersection pattern obtained by cutting through 
a three-dimensional Lorentzian triangulation is a simplicial manifold consisting 
of triangles and rectangles). These concepts can be extended straightforwardly 
to the case where the base manifold B is higher- dimensional. Because of the 
Lorentzian structure of the triangulation, one of the directions in B will be the 
time direction. For example, in the simplicial analogues of spherically symmetric 
models discussed below, dim[B) = 2, with one time t and one radial direction r. 
Also in this case, there will be a well-defined notion of an /-dimensional fibre 
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F over an arbitrary non- integer base point {t,r). This fibre takes the form of a 
piecewise flat manifold, with generalized flat building blocks. 




Figure 2: An edge tower ct x F of a two-dimensional Lorentzian triangulation can 
be represented by a one-dimensional graph consisting of long and short edges which 
corresponds to a cut through the tower at some time which does not coincide with 
t = 1/2. In the figure, the cut at t < 1/2 is indicated by the dashed line, the short 
edges have length /„ and the long ones Id- Because of our planar representation, the 
horizontal distances in the interior of the strip are again not rendered faithfully. 

An important and useful observation is the fact that the d-dimensional geom- 
etry of a cr-tower, where cr is a 6-simplex in B, can be deduced entirely from the 
geometry of the fibre F over a single interior point of a. A simple example is 
again given by a two-dimensional Lorentzian triangulation. The fibre of constant 
time s over some interior point s of a strip [t, t+1], with t an integer, is a sequence 
of straight edges. In order to be able to distinguish the edges that come from 
cutting up- and down-triangles, we must choose s ^ t+1/2 (a similar exclusion 
of symmetric points applies also in more complicated examples). This will result 
in edges of two different lengths and Z^, as depicted in Fig. |2l and it is obvious 
that we can reconstruct the entire strip (or cr-tower) from the knowledge of the se- 
quence of the two types of edges occurring in the fibre over s. An alternative way 
of keeping track of the two different types of edge that can occur in the tower is to 
colour-code them (Fig. El). Both procedures generalize to more complicated prod- 
uct manifolds, with each a-tower represented uniquely by either a (generalized) 
/-dimensional triangulations or a special type of multicoloured graph (see the 
following sections for further examples, and the appendix for geometric details), 
and the matching conditions for a-towers taken into account appropriately. 
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Figure 3: Alternatively, the same edge tower can be represented by a one-dimensional 
graph consisting of edges of two colours, which can be thought of as being induced 
from colouring the up- and down-triangles differently and considering a cut at midtime 
t = 1.5. 

2.1 Lorentzian product triangulation with two-dimensional base 

In preparation for later sections, we will now describe the geometry of Lorentzian 
triangulations T with a two-dimensional base space, for d = 2 + l and d = 2 + 2. 
In both cases, the top-dimensional simplices of B are two-dimensional triangles, 
and T consists of triangle towers. By assumption, the base manifold is itself of 
the form of a (H-l)-dimensional Lorentzian product triangulation B = B' x F', of 
the type usually considered in two-dimensional causal dynamical triangulations 
PP and depicted in Fig. ^ above. Its base B' is a one-dimensional triangulation 
consisting of timelike edges a', with associated edge towers a' x F' consisting 
of sequences of up- and down-triangles in the spatial direction. For our general 
discussion, it will not play a role whether the topology of these strips is spatially 
open or closed. Note that because of the physical interpretation of the triangula- 
tions as causal spacetime geometries, there is no symmetry between the time and 
spatial direction of B. Since the number of triangles in an edge tower a' x F' is 
in general a function of time, B cannot be thought of as a product triangulation 
with the spatial direction as its base. 

It follows that our simplicial manifolds have the form of "staggered" product 
triangulations, which can be thought of as fibrations over B or B'. Accordingly, 
their d-simplices can be characterized by how they fit into either of these product 
structures. The information concerning the fibration T = B' x F' is simply the 
number pair [ii,«2] which counts how many vertices of the simplex lie in either 
one of the two adjacent slices of constant time t, as described above. From this 
specification one can uniquely deduce the geometry of the d-simplex, i.e. which 
of its edges are spacelike and which are timelike, and compute all volumes and 
angles of the simplex and its subsimplices. Analogously, the orientation of a d- 
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simplex with respect to the fibration T = B x F is specified by a triple of numbers 
bi) J2; is]; which count how many of its d+1 vertices lie in each one of the three 
vertex towers contained in a x F. 

2.2 The example of 2+1 dimensions 

In order to train our geometric imagination, we start with the simpler case d = 3 
where the fibres F are one-dimensional. From the point of view of the time 
fibration T = B' x F', the three-simplices or tetahedra come in three different 
time orientations, [3, 1], [2,2] and [1,3]. From the point of view of the fibration 
T = B X F, there are three possibilities how a tetrahedron can appear in a 
tower above a given base triangle a, which in an obvious notation are labelled by 
[2, 1, 1], [1, 2, 1] and [1, 1, 2]. Consequently, one can visualize the tower over a as 
a prism with triangular base, which itself is of the form of a linear sequence of 
tetrahedra of these three types. Assigning three different colours to the differently 
oriented tetrahedra, we note that the geometry of the triangulated prism can be 
encoded in a one-dimensional graph made of edges of the three colours (Fig. 
namely, the graph corresponding to the fibre over any interior point of a. This is 
a three-dimensional analogue of the situation depicted in Fig. El 




Figure 4: A triangle tower of a (2+l)-dimensional product triangulation can be repre- 
sented by a one-dimensional graph consisting of edges of three colours, corresponding 
to a fibre over an interior point of the triangle. The three colours correspond to the 
different orientations with which a tetrahedron can appear in the prism. The figure 
depicts a sequence of tetrahedra [2, 1, 1], [1, 2, 1], [2, 1, 1], [1, 1, 2], etc., where the labels 
of the three vertex towers are as indicated. 

In order to formulate the matching conditions between neighbouring triangle 
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simplex type \ in graph with colours 


111212, 




k 


edge of colour ii 








vertex 


[111 5 5 -'-'■s] 


[1*1? 1*2? 


[1*1 ) , 0,3] 



Table 1: From this table, one reads off the unique building block [jii, jij, jia] in 
the triangle tower of a (2+l)-dimensional product triangulation corresponding to 
a simplex of a one-dimensional coloured graph (left column). The building block 
depends on the fibre type, as expressed by the number of colours occurring in the 
coloured graph to which the simplex belongs (top row). The three ik are colour 
labels. 

towers in this representation, one needs to translate the geometry of the edge 
tower common to the two triangle towers into a one-dimensional graph. One 
simply considers the fibre over an interior point of the edge in B in question. 
It is clear that this gives rise to a two-coloured graph, because it represents a 
cut through a two-dimensional triangulated strip. Denoting by ai and (J2 the 
two triangles in B adjacent to the edge, such a two-coloured graph is obtained by 
moving an interior point in one of the cxi's toward the shared edge. In the process, 
one type of coloured edge disappears from the fibre above the point, because the 
corresponding type of tetrahedron shares only an edge with the common edge 
tower, and not a triangle. The matching condition can therefore be formulated 
as a condition on the linear graphs characterizing o"i x F and a2 x F with one 
colour each deleted. 

We close this subsection with some further geometric observations, which are 
straightforward in a three-dimensional context, but generalize to the more difficult 
four- dimensional case. First, one may also associate a linear graph to a vertex 
tower, which can only ever consist of edges of a single colour. Second, to this 
graph and any other graph corresponding to a fibre over a point of a the Euler 
relation 

No = N,+p (1) 

must apply, relating the number A^^q of vertices in the graph to the number A^i 
of edges, irrespective of their colour. The variable p is either or 1, depending 
on whether the graph topology is or [0, 1]. Tabled lists which building block 
in the triangle tower corresponds to which element (vertex or coloured edge) of 
a one- dimensional coloured graph, depending on whether the graph comes from 
a fibre over an interior point of a (and thus has three colours ii, 12-, is), from a 
fibre over an edge (excluding the vertices) of cr (with two colours, e.g. ii, 12), or 
from a fibre over one of the three vertices (with a single colour) , e.g. i). For an 
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interior point of a, relation (Q) implies 

Niu = N211 + Ni2i + + P, (2) 

where the subscripts of the counting variables Nj^j^j^ refer to the characterization 
of tetrahedra in the triangle tower according to the numbers [ji,j2,i3] of their 
vertices that lie in the vertex towers 1, 2 and 3. Subsimplices of a given tetra- 
hedron inherit a triplet [ji , j2 , Js] in an obvious way, with the three ji adding up 
to the dimension of the subsimplex plus one. This accounts for the appearance 
of the counting number of triangles on the left-hand side of equation Lastly, 
note that for the set of graphs associated with a given base triangle, edges of 
a particular colour occur with the same multiplicity in all graphs if the colour 
occurs in the graph at all. This translates into the relation 

iV2ii = iV2io = N^m = N200 (3) 

and permutations thereof. It follows that all numbers Nj-^j^j.^ in one triangle 
tower are determined by A^2ii) ^121 and A'^ii2- Moreover, these three numbers are 
independent, since each type of tetrahedron can be inserted at any position of a 
triangular prism. 

2.3 The example of 2-\-2 dimensions 

In comparison with the previous subsection, we will use the same simplicial base 
space B, but increase the dimension of the fibre F from one to two. Let us call 
the spatial direction of the base space the "radial direction" (as will indeed be the 
case in later sections). The four possible types of four-simplices from the point of 
view of the time fibration are [4, 1], [3,2], [2,3] and [1,4]. These building blocks 
can still have different orientations within a given triangle tower. Without loss 
of generality, let us assume that the three vertices of the base triangle a have 
coordinates (ti,ri), (ti,r2) and (t2,rs). For a [4, 1] -building block, say, there 
are now three possibilities how the four vertices in the vertex tower over the 
time ti (with respect to the fibration B' x F') can be distributed over the two 
vertex towers (ti,ri) and (ti,'r2) (with respect to the fibration B x F), labelled 
by [3,1,1], [1,3,1] and [2,2,1]. Similarly, a [3, 2] -building block can occur with 
two orientations, [2,1,2] and [1,2,2], but there is only a single way, [1,1,3], to 
orient a [2, 3]-building block. 

The question now arises of how the geometry of a triangle tower can be cap- 
tured by looking at the two-dimensional i?-fibres over (interior) points of its base 
triangle a, which are cuts through T of constant time and radius. The top- 
dimensional building blocks of these surfaces are triangles (coming from [3, 1, 1] 
and its permutations) and rectangles (from [2,2,1] and its permutations). Like 
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sim])lex ty])e 


D' X F' ty])e 


fil")re F 


fil)ro F 


dual gra])li 


[3.1.1] 


[1-1] 


tetraliedroii 


triangle rrr 


■3-\'al. A'ortox r 


[1,3, 1] 


[4,1] 


tetrahedron 


triangle ggg 


3-val. vertex g 


[1,1,3] 


[2,3] 


prism 


triangle bbb 


3-val. vertex b 


[2, 2, 1] 


[4, 1] 


tetrahedron 


rectangle grgr 


4-val. vertex gr 


[2,1,2] 


[3,2] 


prism 


j_ 111 

rectangle rbrb 


4-val. vertex rb 


[1,2,2] 


ro ol 
[3,2] 


prism 


rectangle gbgb 


4-val. vertex gb 


ro 1 1 1 
[2, l,ljr 


ro 1 1 
[3,1] 


triangle r 


edge r 


edge r 


[l,2,l]g 


[3,1] 


triangle g 


edge g 


edge g 


[1,1, 2]b 


[2,2] 


rectangle b 


edge b 


edge b 


[1,1,1] 


[2,1] 


edge 


vertex 


polygon 



Table 2: Four-simplices in the triangle tower over a base triangle with vertex 
coordinates (ti,ri), (ti,r2) and {t2,r3), how they appear in the product triangu- 
lation B' x F' with time base, and in the fibres F' and F, together with a dual 
representation of the latter. The three colours are indicated by r, g and b. 

in previous examples, we can now colour-code subsets of edges in F that will 
always appear with a common length, independent of the base point in B of the 
fibre. Again they fall into three sets, which we will associate with the colours 
red, green and blue. The triangular building blocks in F are monochrome, while 
the rectangles are all bi-coloured, with opposite (and parallel) sides of identical 
colour. The fibre F therefore takes the form of a piecewise fiat manifold with 
triangles and rectangles which are glued together along edges of identical colour. 

Alternatively, it is sometimes convenient to use the graph dual to this (general- 
ized) triangulation. Instead of triangles we then have dual monochrome trivalent 
vertices, and instead of the rectangles dual bi-coloured four-valent vertices, con- 
sisting of pairs of mutually crossing edges of different colour. Since the subgraphs 
of a single colour close on themselves, the dual graph has the form of a superpo- 
sition of three trivalent monochrome planar graphs. Similar to what happened in 
the (2 + l)-dimensional example, if the base point of the fibre is an interior point 
of an edge of cr, (dual) edges of one of the colours disappear, and if the base point 
coincides with a vertex of a, there are only (dual) edges of a single colour left. 

Summarizing our findings, we can say that the geometries of triangle, edge 
and vertex towers of a (2 + 2)-dimcnsional product triangulation are uniquely 
characterized by superpositions of three, two and a single monochrome trivalent 
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building block \ in dual graph 




iii2 


^1 


3-val. vertex ii 


[3. 1 . 1.1 








4-val. vertex iii2 




[2ii , 2j2 , Ojg 




edge ii 


[2ji , li2 5 Ija] 


[2ii , li2 5 Ojg 




polygon 


[lii ) li2) 1*3] 


[111 ) li2) 0^3] 







Table 3: From this table, one reads off the unique building block [ji-^^, ji2, jis] in 
the triangle tower of a (2 + 2)-dimensional product triangulation corresponding 
to a building block of a dual planar coloured graph (left column). The building 
block [jii, ji2, jia] depends on the fibre type, as expressed by the number of colours 
(three, two or one) occurring in the coloured graph to which the simplex belongs 
(top row). The three ik are colour labels. 

planar graph, respectively.^ Table |21 summarizes the simplex types that can occur 
in a triangular tower, and how they appear with respect to the various fibrations. 

With this characterization in hand, we can now apply the Euler relations for 
planar graphs, 

- m +m = x 

3N^^'^* + 4A/'^^^* = 2N*, (4) 

relating the numbers A^^ of building blocks of dimension d contained in the dual 
graphs associated with a triangle tower. In (^, NQ=Nlf^* + N^'^^* is the sum of 
dual three- and four-valent vertices, N* and A^2 ^^e numbers of dual edges 
and faces, regardless of their colour, and x is the Euler number of the fibre. 
Because there are different graphs associated with base points in the interior of 
the triangle a, and its edges and vertices, (jH) amounts to a total of 14 equations 
for the counting numbers Nj^j^ja the various simplex types of a given triangle 
tower. They can be used to express and functions of n'^'^* and N^f^* , 

say. Applying Q to the graph representing the full triangle tower and making 
use of the "translation table" , Table 01 we obtain 

A^311+A^131 + A^113 + A^221 + A^212 + A^122-A^211-A^121-A^112 + A^lll = X, 
3{N3ll + Ni31 + Nu3)+^{N221 + N212 + N122) = 2{N2U + Ni21 + Nu2). (5) 

^This picture is reminiscent of three-dimensional causal dynamical triangulations, where the 
three-dimensional 'sandwich' geometry of a discrete time step At = 1 can be represented by a 
dual graph which is a superposition of two monochrome trivalent planar graphs |2()| . 
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Furthermore, the numbers of three- and four-valent vertices of the same colour(s) 
are equal for each of the graphs in which they appear, that is, 

^311 = ^310 = ^301 = ^300; 

N221 = N220, (6) 

as well as permutations hereof. In conjunction with eqs. (^, they imply that each 
^hhi-i is determined by the six numbers iVsn, iVisi, A''ii3, N221, N212 and A^i22- 
Moreover, within each triangle tower, these numbers are independent, because 
one can alter them separately by applying local changes in the geometry. 

3. Extrinsic curvatures and light expansion rates 

The main aim of this section is the derivation of an expression for the light 
expansion rate of a spacelike (rf — 2)-dimensional surface 5* of constant time 
and constant radius in a simplicial product manifold of dimension four, in order 
to formulate necessary criteria for the presence of black holes in the quantum 
theory. We start by reminding the reader of the geometric meaning of trapped 
surfaces. For reasons of completeness, we then review the construction of the 
expansion rate in terms of the extrinsic curvatures of S and a suitably chosen 
spacelike hypersurface S for the case of a smooth c?-dimensional manifold (M, Qab), 
following the treatment in [21]. We then translate this to simplicial manifolds, 
starting with a construction of the extrinsic curvatures. The expansion rate is 
first computed for a three-dimensional triangulation, because the geometry of 
the situation is closely analogous to that in four dimensions, but much easier 
to visualize. We discuss the role of different types of simplicial building blocks 
and their effect in terms of the focussing and defocussing of light rays. As one 
would expect, the simplicial expressions for the expansion rates are subject to 
discretization ambiguities. We find that for specific choices of how the piecewise 
flat surface S traverses the top-dimensional building blocks of the triangulation, 
the formulas for the expansion rates in both three and four dimensions take on a 
particularly simple form which will be useful in numerical simulations. 

3.1 Trapped surfaces 

A central question in our attempt to set up a path integral for black-hole geome- 
tries is how one may recognize the existence of a black hole region or a horizon 
in the quantum geometry. One difficulty is that the usual definition of a black 
hole region is both a classical concept and highly non-local (for details, see j22I 
and references therein). To find the black hole region one has to know the entire 
spacetime manifold A4 = {M,gab) and to find the entire causal past J~{'^~^) of 
future null infinity 5^. The spacetime Ai is then said to contain a black hole 
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region if Ai does not coincide with the causal past of future null infinity. What 
motivates this definition is the fact that "nothing can escape from a black hole" , 
and, in particular, nothing can escape to future null infinity. The event horizon 
is defined as the boundary of the black hole region, which in turn is defined as 
the complement of J^{Z^) in Ai. 

This definition is not suited for the formulation of a nonperturbative path 
integral in terms of dynamical triangulations, because there one can only work 
with spacetimes of finite volume. A more local criterion for a horizon or a black 
hole region is the existence of a so-called trapped surface. This also appears 
in other approaches to quantum geometry (see and references therein), in 
classical general relativity in order to define "(nearly) isolated" horizons j21], as 
well as in "horizon finders" in numerical relativity (see, for example, [2S])- 




Figure 5: Inward- and outward-pointing light rays (solid and dashed arrows) outside a 
black hole region (left) and at an apparent horizon (right). Here we depict the situation 
in a (2-|-l)-dimensional spacetime, where the surface S (lower circle) is one-dimensional. 

A trapped surface 5* is a compact spacelike surface of codimension two, with 
the following property. Consider future-directed null geodesies (i.e. light rays), 
which start off orthogonally from S. This can happen in two directions, called 
"inward-pointing" and "outward-pointing" , see Fig. |S1 (The distinction is natural 
since S is assumed to be closed, for instance, a sphere.) Light rays which point 
outward, say, can either diverge from or converge toward each other locally. This 
effect is measured by the outward-pointing expansion rate which can be either 
positive (for divergence), negative (for convergence) or vanishing. For a sphere 
in fiat Minkowski space the expansion rate if+ is everywhere positive and the 
analogous expansion rate H_ for inward-pointing light rays negative. However, 
for a sphere inside a black hole, even the outward-pointing light rays are bent 
inwards to such an extent that the expansion rate if+ becomes negative too. 

Consequently, one defines a trapped surface as a compact spacelike surface of 
codimension 2, whose expansion rates if+ and if_ are both negative. If "neg- 
ative" is replaced by "non-positive" one speaks of a marginally trapped surface. 
The outermost marginally trapped surface is called the apparent horizon (assum- 
ing certain technical conditions, which are unimportant here). Its geometry is 
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depicted in Fig. One can show that its expansion rate vanishes every- 
where. 

For a Schwarzschild solution the apparent horizon coincides with the event 
horizon. Moreover, a marginally trapped surface is always contained in a black 
hole region, as defined above (see [221 ^^"^ references therein for a precise for- 
mulation of the theorem). This holds also for spacetimes which do not fulfil the 
Einstein equations, but instead the condition Rabk"'k^ > for the Ricci tensor 
and all null vectors k"'. The latter is satisfied if the Einstein equations hold and 
if the matter satisfies the strong or weak energy condition. 

These classical considerations motivate the replacement of the global charac- 
terization of black holes in terms of event horizons by the more local criterion 
that a trapped surface exist. The condition Rabk"'k^ > will not in general be 
satisfied by the spacetime geometries contributing to the nonperturbative path 
integral, because the individual path-integral histories are arbitrarily far away 
from classical solutions. However, one would expect to recover such a property 
at sufficiently large scales in an appropriate continuum limit of the theory. We 
will in the following concentrate on the condition = as an indicator for the 
presence of black holes. 

3.2 Continuum treatment 

The light expansion rate H describes the expansion of a family of lightlike geodesies 
starting orthogonally from the (rf — 2)-dimensional surface S. It is defined by 

H = VaV". (7) 

where V is the Levi-Civita connection associated to the spacetime metric gab 
and p° is the tangent vector field to a future-pointing null geodesic congruence 
starting off orthogonally from the surface S. In other words, p"" must satisfy 

PaP" = 0, rUabP'is = 0, (8) 

and 

p'^VaP' = 0, (9) 

where rriab is the metric induced on 5* by the spacetime metric gab- Embedding 
S into a spatial (constant-time) hypersurface S with future-pointing unit normal 
vector field n"' and calling ±g" the (out- and inward-pointing) unit normals to S 
tangential to E, the ((i— 2)-dimensional metric is given by 

rriab = gab + naTib - qaQb- (10) 
For later use, we also introduce the notation 

hab = gab + naUb = rUab + qaQb (11) 
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for the induced metric on E. Note that niab is independent of the choice of S as 
long as the latter contains S. The general solution to eqs. (jHI) on S" can then be 
written as 

±pf5 = c±K±g«), (12) 

with coefficients c± which are positive functions on S. As long as the functions 
c± are kept arbitrary, the general solution (fT^ does not depend on S in the sense 
that if one starts with another hypersurface S' and corresponding normal vectors 
n'"' and g'", one can always find functions c'_j_ such that 

15 = c±{n^ ± Q") = c'^in"" ± q"") = ^p']s- (13) 

Using the decomposition (fTUI) . we can compute the expansion rates H± for the 
in- and outward pointing light rays, 

= m'^'Va^Pb + qy'Va^Pb - n'^u'Va^Pb 

= m'^'Va^Pb - {{Ci)'' V T q") {c±r' [Vv/pJ ± q^c^y' [^p'Va^Pb] 
= m'^'Va^Pb. (14) 

Going from the second to the third line, we have used eq. (jl2|) . The terms in 
square brackets vanish by virtue of relations (jHl) and Noting that the covariant 
derivative is projected onto the surface 5" in the last line of (HH), we can simplify 
the expansion rate further by using again the expression (fT^ for the vector fields 
on S, 

H± = m'^''VaC±K±gfe) 

= c^m^^VaUb ± c±m"''Vagfe + [m''\nb ± qb)] VaC± 

= c± {-m'^'Kab T k) . (15) 

The term in square brackets vanishes by definition of pb, and in the third line we 
have defined the extrinsic curvature 

Kab = -KV.Ub (16) 

of S in M, and the extrinsic curvature 

kab = -m^X^aqb (17) 

of S in E. 

The last line of (jl5|) shows that the signs of the expansion rates H± on S 
(which determine whether or not S* is a trapped surface) are independent of the 
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prefactors c± introduced in eq. (fT^ which are always positive. Likewise, the 
condition if+ = for an apparent horizon does not depend on the choice of the 
constant-time surface S and therefore on c+. Our reason for keeping track of 
the prefactors c± exphcitly is the fact that in the quantum theory it may be 
convenient to monitor integrated expansion rates rather than local ones. If such 
an integration were performed over a surface S* at a given time and radius and 
if the underlying geometry were exactly spherically symmetric, the expansion 
rate would be constant on S and the S-dependence of the integral would simply 
amount to an overall factor c±. However, in situations without exact spherical 
symmetry, either in a smooth or a simplicial setting, it does matter a priori that 
the absolute magnitude of H± depends on c±. If we define an "averagely (mini- 
mally) trapped surface" 5* in a generic geometry and for a specific hypermanifold 
Eby 

S)+{S):= [ H+VdetmdS = 0, (18) 
Js 

where V det m dS is the invariant volume element on S, the same condition will 
in general not hold for a different choice S' of the hypermanifold containing S, 
as was pointed out in [2^1 ■ The average expansion rate S^'^ may become positive, 
say, because the prefactor c'^ which appears when one expresses = -|- q"- in 
terms of the primed normals, 

= c'^{n"' + q"'), (19) 

may be large in a region of positive (defined with respect to S) and small in a 
region of negative whereas the contributions from the two regions cancelled 
each other with respect to the unprimed spatial slice S. 

It follows that some caution has to be applied when using the vanishing of 
averaged expansion rates as an indicator for the presence of black holes. This ap- 
plies in particular to the case of simplicial manifolds which by their very nature 
can never be exactly spherically symmetric. Nevertheless, there are important 
situations - including the one considered here - where the vanishing of the inte- 
grated expansion rate, ()18|) . is the relevant criterion for indicating the presence 
of an apparent horizon in the quantum theory. Firstly, we may use the preferred 
time foliation shared by all Lorentzian triangulated geometries in a quantum su- 
perposition to define the null geodesic congruence in terms of eq. (jl2j) with c± = 1 
and thus arrive at unique values for H±. Second, as explained in the Introduc- 
tion, we are interested in studying the path integral for approximately spherically 
symmetric geometries, in which case the variations across S of the expansion rate 
will not be large. Moreover, the surfaces S will be approximately spheri- 
cally symmetric, and therefore the factors c'_j_ connecting the expansion rates of 
two different (approximately spherically symmetric) surfaces S and S' will be 
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approximately constant, so that the case described above will generically not 
occur. 

One can also formulate (local or integrated) criteria for the presence of trapped 
surfaces which depend quadratically on the expansion rates H±. For example, a 
potential advantage of using the product is the fact that it is independent 

of the choice of the constant-time surface E j^E]- A generic problem with such 
expressions in a simplicial approach is that they contain products of delta func- 
tions at the same point, and any regularization procedure is subject to a high 
degree of discretization ambiguity. Also, one generically cannot avoid an explicit 
dependence of the formulas on the local geometry of the triangulation, as opposed 
to mere counting of simplex types, which can be achieved for the linear expansion 
rates (c.f. Sees. 3.4 and 3.5 below). For these reasons, we will not pursue this 
possibility presently. 

3.3 The case of piece wise linear manifolds 

In this section we will develop an expression for the expansion rate for the case of 
a piecewise fiat manifold, and will assume that both the constant-time surface S 
and the surface S are fibres with respect to the two product structures involved. 
The second of these assumption we will relax later on in Sees. 3.4 and 3.5, for 
reasons explained there. In order to apply formula (jl5j) to the case of a piecewise 
linear manifold, we have to define an expression for the extrinsic curvature of 
a (piecewise linear) hypersurface in such a manifold. We begin by considering 
the extrinsic curvature Kab of a (rf — l)-dimensional constant-time hypersurface 
Sj (see the appendix for a definition) in a c?- dimensional triangulation. Since 
the geometry of each ci-simplex is fiat and Minkowskian, a constant-time surface 
inside it is just a linear spacelike hypersurface (with boundaries) of this Minkowski 
space, and its extrinsic curvature vanishes. We will arrive at the same result (jHUj) 
as reference ^H] • However, our prefactors will differ from the ones obtained there, 
because we use different coordinates, which also appear as argument of the delta- 
function. 

A key observation is that one can always imbed two adjacent d-dimensional 
simplices af and cxg isometrically into a common (i- dimensional Minkowski space. 
It follows that the geometry across the (rf — l)-dimensional boundary simplex 
^ira separating the two adjacent simplices is fiat. (This no longer holds for the 
(d— 2)-dimensional subsimplices, on which the intrinsic curvature is concentrated.) 
In the common coordinate system, light rays crossing the (d — l)-dimensional 
timelike boundary simplex appear as straight lines, but a constant-time surface 
will in general be seen to have a "kink" there. This implies that its extrinsic 
curvature is concentrated on the (ci — 2)-dimensional intersections fl af^l of 
the constant-time surface with the (c?— l)-dimensional timelike subsimplices. 
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The relevant geometry is therefore that of two hnear spacehke hyperplanes in 
Minkowski space meeting in a kink (a (d— 2)-dimensional hnear submanifold) with 
a certain relative angle. The geometry is homogeneous in the directions along 
the kink. In the triangulation, we can ignore these directions as long as we stay 
away from the (ci — 2)-dimensional boundary simplices of the (ci — l)-dimensional 
timelike subsimplex. 




Figure 6: The geometry of two neighbouring d-simphces, seen in the two-dimensional 
plane perpendicular to the kink formed by the intersection of a constant-time surface 
(dotted lines) and their common (d— 1) -boundary simplex. The common Minkowski 
coordinate frame is spanned by the unit vectors eg and ei, with corresponding coor- 
dinates t and x^. The origin of the spatial coordinate is taken to coincide with the 
location of the boundary simplex. Note that the normal vectors are perpendicular in a 
Lorentzian sense to the constant-time surface. 

To analyze the geometry of the two-dimensional plane spanned by the two 
vectors normal to the kink, it is convenient to use the basis vectors (eo)° and 
(ei)'' orthogonal to ("1(7^^2 (see Fig. IH)). By convention, (eo)° is future-directed 
and parallel to (jf^2 i^i)"" points from af to erf. Furthermore, we introduce 
corresponding Minkowski space coordinates which measure the proper 

distances along the integral curves of (eo)'^ and (ei)'* respectively. We define 
to be zero on crf|^2 . 

For the computation of the extrinsic curvature tensor we need the unit normals 
n%^ and n%. to the two pieces Tii^^ and Ti?^ of the constant-time surface. Because 
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of their timelike nature, they can be written as 

n^,) = cosh Pi (eo)" + sinhp^ (el)^ (20) 

where pi,i = 1,2 is the hyperbohc angle between n(j) and eg. Similarly, the 

(spacelike) unit tangent vectors to sj"^^ and sl^"* (in the span of (bq)" and (ei)") 
are 

4) = sinh Pi (eo)" + cosh pi (ei)^ (21) 

According to eq. (|16|). the extrinsic curvature is given by the projection onto 
of the covariant derivative of the normal vector field, which in our Minkowski 
coordinates reduces to the usual coordinate derivative. In order to deal with the 
discontinuity of the normal vector field of at the kink, we regularize it with the 
help of a family of smooth functions 5^ which converge to the delta function as 
£ — > 0. The angle p between the time axis and the normal is then approximated 
by 

p{x') = p, + Ap J 6e{x') dx'\ (22) 

where Ap = p2 — pi- Projecting —'Van''{x^) onto the hypersurface with the 
projection operator h'^ := 7]^ + nan^ (where r/ac is the Minkowski metric), we 
obtain the extrinsic curvature 

K^\x^) = -{r]^^ + n''{x^)n%x^))Vc{cosh p{x^) {eof + sinh p{x^) (eif) 
= -{ri"' + n''{x')n%x^))5eix') Ap (ei)e s\x^) 

= -6,{x^) Ap coshp(x^) s^ix^) s\x^), (23) 

where we have used that VaX^ = (ei)^. 

For the calculation of the extrinsic curvature k"''' of the (rf — 2)-dimensional 
surface S in the surface we proceed completely analogously. The constant-time 
surface is a piecewise linear manifold, albeit a generalization of a simplicial 
manifold with more general building blocks than simplices. Again, we can now 
embed any two adjacent ((i— l)-dimensional building blocks of the Sf-triangulation 
in a (d— l)-dimensional fiat euclidean space. We choose the surface 5 as a linear 
(o? — 2)-surface inside each building block, possibly with a kink when one crosses 
from one building block to the other. 

To analyze the geometry of S, consider the same two (i-simplices as before 
and add a third basis vector e2 which is defined to be orthogonal to eo, ei and 
to the intersection S fl ai^2- d = 3 the intersection is just a point, and for 
(i = 4 the geometry is homogeneous along the intersection. In the latter case, we 
complement the basis with a fourth normalized vector 63 orthogonal to the other 
three and thus tangential to 5* fl o-f^^g . 
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x^ = 



Figure 7: The geometry inside a constant-time hypersurfacc S^. Tlic figure illustrates 
various geometric quantities wliich appear in the calculation of the extrinsic curvature 
of the surface S (dotted line) in S^. The surface S has been chosen straight inside 
either of the two neighbouring c/-simplices, so that the only nontrivial contribution to 
the extrinsic curvature comes from the kink of S on the boundary between the two 
simplices. The vectors S(i) and S(2) may differ in the direction of eo which is suppressed 



where the angle ■0^ is the angle between the normal and the unit vector 62- 
Repeating our previous construction, we define an interpolating angle ip by 



with = -02 — ^1- To determine the extrinsic curvature of S, we begin by 
calculating 



here. 



The (outward-pointing) unit normals q'^^^ to S'^ :— S (1 af in and those 
tangential vectors u^^-^ to S'^ which lie in the plane spanned by 62 and may be 
written as 





(25) 




(26) 
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where denotes the induced covariant derivative in the submanifold E^. This 
leads to the expression 

= -6,{x^) coshp(xi) cos^/'(xi) u^{x^) u^x^). (27) 

for the extrinsic curvature of S* in S^. Given the extrinsic curvatures of and 
S, eqs. and (j27j), it is now straightforward to apply the continuum formula 
(fT^ (with c+ = 1) and calculate the expansion rates of outward-pointing light 

rays as 

H+{x') = -k-m^'Kab 

= 6s{x'^){^ coshp(x^) cosipi^x^) + Ap coshp(a;"^) cos^ ■?/'(a;^)).(28) 

For the computation of the integrated expansion rate, it is convenient to 
introduce the variable p measuring the proper length along the orbits of m". From 



u"- = - sin V^(e2)" + cos V'(sinhp(eo)" + cosh p(ei)'') (29) 

we have dx^ = cosh p cos ip dp, and it follows that 

H+ (p) = 6e{p){^ + Ap cos ij{p)) (30) 

"'^Uip)i/klj + Apcos^lJm) (31) 



where ipm '■= |('^i + "^2) and we have assumed the functions Se to be symmetric. 
For the case of a four-dimensional triangulation, by introducing a second (proper- 
length) coordinate x^ along (63)° the invariant volume element V det m dS of S 
assumes the simple form dp dx^ . We can now easily integrate the expansion rate 
over a neighbourhood in S enclosing the set S fl aff^2 (that is, the support of the 
expansion rate in S), resulting in 

[ H+dp dx^ = Vol((Tfn2 n S) {At/j + Ap cos ^„) . (32) 

inbh. 

(For a three-dimensional triangulation, d = 3, we set Vol(crfpi2^ fl S') = 1). In the 
integral over the whole surface S, the contributions from the (d — 3)-dimensional 
building blocks simply add up to give the integral iD+(5') of the expansion rate 
over the surface S, 

9)+{S) = / H+VdetmdS = ^ Vol(^d-in5)(A/' + A^cos V'm)(^d-in5)- (33) 

When evaluating an expression like (j33p for a dynamical triangulation it is 
inconvenient that the angles pi and ipi defined by eqs. (pUj) and depend on 
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the indexing of the d-simphces cJi and a2- This motivates the definition of new 
angles pi = (— and ipi = {—lyipi, in terms of which we have 

cosh Pi = -{eo)""n\i)r]ab, sinhp^ = (§1)^.) nj.)?]^^, Ap = pa + Pi (34) 
cos = (62)" gfi) r]ab, sin ijji = 5^.) gj*-) r]ab, Ai/' = T/'a + V'l, (35) 

where (ei)^-j = (— l)*(ei)" and s"-^ = (— l)*s"j-| are pointing toward the simplex 
erf for both 2 = 1 and z = 2. Furthermore, cosipm may be written as 

cos V„ = cos(i(Vii - V^a)) = cos(i(-^i + ^2)), (36) 

so that all quantities Ap, cos ipm appearing in the integrated expansion rate 
(j33|) are independent of the index i in af. 

3.4 The case of dynamical triangulations: 2+1 dimensions 

In order to get an idea of what is involved in applying the framework of the 
previous section to the particular case of causal dynamical triangulations, we 
will first consider the simpler case of 2+1 dimensions. The triangulations take 
the form of fibrations over a two-dimensional base manifold, parametrized by 
a time and a radial coordinate. Our aim will be to derive a formula for the 
(integrated) expansion rate which is operationally simple to evaluate. In 2+1 
dimensions the surface S on which we want to define the expansion rate if+ is one- 
dimensional and if+ is concentrated on the zero-dimensional vertices v = (1 S 
of the piecewise straight surface S. Following standard convention, we set the 
volume of these vertices to 1. Consider now a fibre S' in a [Ir.m, ^T,out, 1t+i]- 
triangle tower. ^ Applying formula (j33j) for the integrated expansion rate gives 

9j4S) = -i{S)+J2Mv)cos^Uv), (37) 

vCS 

where 

e(5) = |(iV2n-iVi2i) (38) 

is the integrated extrinsic curvature of S in Sj. The values of the angles A^(f ) = 
Pi{v) +p2{v) and cos{iprn){v) = cos(|('?/'i(f ) —ip2{v))) contributing to eq. (jHTf) can 
be read off the Tables 0] and |S1 The term (jHHj) is just a counting term which only 
depends on the total number of simplices of a certain type. This is not true for the 
remaining sum over vertices v in eq. (j37|) . To evaluate it, we need to know not only 

^Here T is an integer-valued time, and the lower-case t will henceforth denote a time value 
between T and T + 1. The subscripts in and out refer to vertices with radial coordinates ri„ 
and rout, with < rout- 
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simplex 


boundary simplex 


P 


[3, 1] 


[2,1] 


— arsinh^== 


[1,3] 


[1,2] 


arsinh^== 


[2,2] 


[2,1] 




[2,2] 


[1,2] 





Table 4: The angles p per simplex type and per boundary simplex type contributing 
to the expansion rate in 2+1 dimensions. As usual in causal dynamical triangulations, 
we fix all spacelike edges to have a squared length l1 = and all timelike edges to 
have l1 = —aa'^, where a is the so-called lattice spacing that will eventually be taken 
to zero. The number a parametrizes the (fixed) ratio between the length units in time 
and space directions. 



[^m,T, jo«t,r, ^Titl] 




[2, 1, 1] 


n 
6 


[1,2,1] 


TT 

6 


[1,1,2] 






Table 5: The angles ip per simplex type in a triangle tower of type [lT,m, ^T,out, 1t±i] 
contributing to the expansion rate in 2+1 dimensions. 

the numbers of the various simplex types but also how they are glued together. 
In other words, there exist pairs of triangulations with identical numbers N^y^ 
which nevertheless have different (integrated) expansion rates. This dependence 
on the local gluing information is caused by the factor cos-^m in (|Tfj). and can be 
traced back to our particular choice of the surface S. In the previous subsection, 
we assumed S to be straight inside each (constant-t hypersurface of a) three- 
simplex, as illustrated in Fig. [7| This leads to a coincidence of angles p and at 
the intersection point S fl af^2 between two three-simplices af and a^. We will 
adopt a different choice for the surfaces S which will simplify the functional form 
of the integrated expansion rate. We expect that this will speed up an eventual 
numerical implementation and also simplify any analytic considerations of the 
evaluation of the expectation value of the expansion rate in the path integral. 
The precise choice of S* is a typical discretization ambiguity, which should not 
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[2,1,1] [1,2,1] [1,1,2] 

Figure 8: The "building blocks" which constitute our surfaces S of constant radius 
(dashed lines) inside the hypersurfaces in 2+1 dimensions. 

make any difference when we perform the continuum hmit7 In the absence of 
an argument for which discretized choice of a "surface of constant radius" is 
more natural, we will adopt the simplest prescription, both here and in the four- 
dimensional case discussed in the next subsection. 

Our alternative hypersurfaces S lie still within Sj, and are piecewise straight, 
but we arrange their straight segments to be dual to the edges of the original 
triangulation in E^, as depicted in Fig. |H1 so that the "kinks" of S no longer 
coincide with the points S fl (jf^2- Since the angles Aip are still the same as for 
the original surface S, the first term in the expansion rate ^7\\ is unaffected, 
t{S) = t{S). However, since cosipm is now always equal to 1, we obtain for the 
integrated expansion rate of 5* 

MS) = -I (iV2ii - iVi2i) - {2arsinh^=i== {N^n + N^i) } + 
o V3 + 12a; 

2 arsinh^^= Nn2 (39) 
VI + 4a 

for a triangle tower of type [lT,m, It.om*, Ir+i]- For a triangle tower of type 
[It, lT+i,m, ^T+i,out] the first term, i.e. the extrinsic curvature of S, is unchanged 
(apart from relabelling), whereas the terms coming from the extrinsic curvature 
of St change sign, leading to 

MS) = -I {Nui - iVii2) + {2arsinh^=i== (iVi2i + Nn2) }- 
o ^/6 + Via 

2arsinh^^= Ar2ii. (40) 

One thing to note about eqs. (jH^ and (pn|) is that the two terms in curly brackets 
have merely boundary character in the sense that they cancel when one considers 

''Obviously such an assertion must eventually be verified. 
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Figure 9: The effect of different tiiree-simplices in a [lT,in, ^T,out, lT+i]-triangle tower 
on future-directed light rays. We only show the angular behaviour, which is relevant for 
calculating the expansion rate. Tetrahedra of type [3,1] focus light rays in the angular 
direction (left figure), whereas [2,2]-tetrahedra defocus them (right figure). 

the expansion rate for more than one time step. For example, suppose we added 
the integrated expansion rates associated with two base triangles of type [1t-i, 2t] 
and [2t, 1t+i] with a common spacelike edge. Then the contributions from the 
two terms in curly brackets cancel each other, because the numbers of [1t-i, 3^]- 
tetrahedra in the first triangle tower equals the number of [3^, l^+i] -tetrahedra 
in the second. 

A second observation about the counting formulas (jHUI) and (jiUI) is that they 
have a direct geometric interpretation in terms of the focussing of light rays pass- 
ing through the three-dimensional building blocks of the product triangulation. 
Let us call the fibre direction angular and the other spatial direction radial, and 
consider first a [lT,m, ^T,out, lT+i]-triangle tower. It contains both [3,1] -tetrahedra 
(which can be split into the subtypes [2,1,1] and [1,2,1]) and [2, 2] -tetrahedra 
(of subtype [1,1,2]). The first term in (jH^ . which depends on the difference 
(Ai'211 — A''i2i), has a purely spatial origin and is determined by the triangulation 
of the spacelike edge tower. The second term (in curly brackets) in (jH^ gives a 
negative contribution, since a [3,l]-tetrahedron focusses light rays in the angular 
direction (see Fig. El left). It also focusses light rays in the radial direction. At 
any rate, since these effects cancel out for adjacent time slices as we have just 
explained, let us concentrate on the remaining term in (jH^ . As illustrated in Fig. 
13 the [1, 1, 2]-tetrahedra defocus light in the angular direction (thus accounting 
for their positive contribution to the integrated expansion rate on S), and focus it 
in radial direction. In a [It, Ir+i.m, lT+i,oMt]-triangle tower the focussing effects 
are exactly reversed: the [1, 3]-simplices defocus light rays (or timelike normal 
vectors) whereas the [2, 1, l]-tetrahedra focus in angular and defocus in radial 
direction. We conclude that apart from the extrinsic curvature term t of the 
(one-dimensional) surface S the expansion rate is determined by simply counting 
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the numbers of focussing and defocussing [2,2]-tetrahedra, associated with the 
[1^-1, 2r] and [2^, ly+ij-base triangles respectively, and taking their difference. 

3.5 The case of dynamical triangulations: 3+1 = 2+2 dimensions 

We will now generalize the treatment of the previous subsection to the case of 
a four- dimensional Lorentzian (i.e. causal) dynamical triangulation which is a 
product of a two-dimensional triangulated "r-t-plane" and a two-dimensional 
fibre, representing the angular directions. We will follow the same strategy as in 
the three-dimensional case to obtain a simple functional form for the expansion 
rate S)+{S) which does not depend on detailed local gluing data. 

For a (2 + 2)-dimensional triangulation, S is two-dimensional and Sj+{S) has 
distributional support on the one-dimensional edges e = H S in S. As before, 
if one defines the surface S" to be a two-dimensional fibre in a triangle tower, it 
turns out that the expansion rate does not just depend on the numbers of the 
various simplex types in the triangle tower, because of the cos ^/'m-terms in eq. 
fl33|) . We will therefore choose an alternative surface S, still inside the constant- 
time hypersurface S^, but which intersects the two-dimensional building blocks 
contained in orthogonally such that all cos ■j/'^-f actors will be equal to one. 

Note that in the (2 + l)-dimensional case the surface 5* was related to the 
dual triangulation of - the vertices of S were positioned at the barycentres 
of the triangles and rectangles which are the piecewise flat building blocks of 
Sj. (The vertices inside the rectangles did not appear explicitly, because they 
do not contribute to the extrinsic curvature of S.) Also in 2+2 dimensions we 
will construct the surface 5* as part of the dual triangulation of the constant- 
time surface Sj. Recall from Sec. 2.3 that Sj consists of equilateral tetrahedra 
and triangular prisms with two equilateral triangular and three rectangular sides. 
To construct the dual triangulation, one places vertices at the barycentres of 
these tetrahedra and prisms, which one then joins by "dual" edges if they belong 
to neighbouring building blocks. These dual edges will orthogonally cross the 
boundary faces (triangles or rectangles shared by adjacent tetrahedra and/or 
prisms) exactly at their barycentres. The dual fiat two-surfaces spanned between 
pairs of dual edges will therefore also intersect the two-dimensional faces of the 
Sj-triangulation orthogonally. Our surfaces S will by definition be built from 
those dual two-surfaces which extend in the angular directions. 

The extrinsic curvature of a surface S with respect to is now concentrated 
on the dual edges. For some dual edges k vanishes, for example, for the dual 
edges in the [2^, 3T+i]-building blocks in a [2t, lT+i]-triangle tower. We will 
now discuss each case in turn and compute its associated angle contribution Aip. 
Fig. ^1 shows the three-dimensional fiat building blocks that can appear in a 
constant-time hypersurface S^, each with its corresponding "elementary" surface 
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out 




out 




Figure 10: The elementary contributions of various three-dimensional building blocks 
in St to the two-dimensional surface S in a [lT,in, ^T,out, lr+i]-triangle tower. The cases 
not shown, associated with four-simplices of type [1,2,2] and [1,3,1], can be obtained 
by exchanging the role of the "in" and "out" vertices of the cases [2,1,2] and [3,1,1] 
respectively. 

contributing to S. There are essentially four different cases: 

(i) The intersection of a four-simplex of type [2,1,2] or [1,2,2] with Ej is a prism. 
The surface S intersects both triangular faces and two of the rectangular 
faces. The dual edges in S which cross the rectangular faces do not carry any 
extrinsic curvature but the two edges which cross the triangular faces do. 
The angular difference between the two normals involved is A-^ — ±7r/3. 

(ii) The intersection of a four-simplex of type [3,1,1] or [1,3,1] with is a tetra- 
hedron. The surface S intersects three out of the four boundary triangles. 
Three dual edges connect the barycentrc of the tetrahedron with the centres 
of these three triangles. The angular difference for each of the three dual 
edges is again At/j — ±7r/3. 

(iii) The intersection of a four-simplex of type [1,1,3] with is a prism. The 
surface S intersects the three boundary rectangles. Since the intersection 
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simplex 


boundary simplex 


P 


[4,1] 


[3,1] 




[3,2] 


[3,1] 




[3,2] 


[2,2] 





Table 6: The angles p per simplex type and per boundary simplex type contributing 
to the expansion rate in 2+2 dimensions. For the time-reversed simplices we have 
P{[i,k]) = -P{[k,i])- 

surface is a single flat triangle parallel to the triangular boundaries of the 
prism, there are no contributions to the extrinsic curvature. 

(iv) The intersection of a four-simplex of type [2,2,1] with Sf is a tetrahedron. 
The surface S intersects all four boundary triangles. However, since there 
are two types of dual edges which contribute with angles of opposite sign, 
the net contribution to the extrinsic curvature vanishes. 

The extrinsic curvature of St with respect to the four-dimensional triangula- 
tion is concentrated on the two-dimensional faces of the S^-building blocks. The 
angles p needed to calculate this quantity are hsted in Table IHl The surface 
S intersects the two-dimensional faces orthogonally, and the extrinsic curvature 
term K"-^mab appearing in the expansion rate (c.f. eq. (|28|)) therefore has distri- 
butional support on these one-dimensional intersections. Note that they do not 
coincide with the dual edges contained in S, but are positioned transversally to 
them. To determine the integrated extrinsic curvatures K we still need to calcu- 
late the lengths of these intersections, which is straightforward and will lead to 
an explicit t-dependence of the expansion rate. Similarly, to obtain the extrinsic 
curvatures k one has to multiply the angular differences Aip with the lengths of 
the associated dual edges. The final result for the integrated expansion rate for 
a [lT,m, Ir.oMt, lT+i]-triangle tower and a surface S" in a constant-time surface 
is given by 

^+{S) = a({l- t)--^(iVi3i - iVan) + t^(iVi22 - N212) + 
^ 2v6 <J 

|(1 _ t)(-v^) arsinh(^— ^=)(iV3n + iVi3i + \N221)] + 
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(41) 



Similarly, for a [1t-i, It,™., lT,oMt]-triangle tower and a surface 5 in a constant- 
time surface one obtains 

Sj+iS) = a(t^{Ni,s - iVi3i) + (1 - t)^(iV212 - iV22l) - 



\t{-V3) arsinh(^— ^=)(iVi3i + iVna + ^Nu2)} 



(42) 



These formulas have the desired form of depending only on counting-variables 
for the simplices in a given triangle tower. As before, one can achieve a further 
simplification by adding contributions from successive time layers. Consider two 
triangle towers which are joined by a spacelike edge in the base manifold. If one 
chooses the two constant-time surfaces at which the contributions to the extrin- 
sic curvature are evaluated to be at times (T + t) and (T — 1) + (1 — t) = (T — t), 
the terms in the curly brackets in (j4ip and (j42j) will cancel each other because 
the numbers of the various types of [4^, lT+i]-simplices equal those of the cor- 
responding 4T]-simplices. An obvious and symmetric choice that achieves 
this cancellation over a larger number of time steps is therefore to evaluate the 
expansion rate always at half-integer times. 

To summarize, the expansion rate in 2-f-2 dimensions is again a sum of terms 
coming from the extrinsic curvature k of S and of the extrinsic curvature K of 
E^. The former are already determined by the triangulation of the spacelike edge 
tower, namely, the partitioning of the [4, l]-simphces into the subtypes [3,1,1], 
[1,3,1] and [2,2,1], and similarly for the [l,4]-simplices. In analogy with what 
happened in the (2+l)-dimensional case, the contribution from the curvature K is 
determined by the distributions of the simplicial building blocks of type [3, 2] and 
[2,3]. In a [lT,m, It.om*, lT+i]-triangle tower the [1, 1, 3]-simplices defocus light 
rays and in a [1t-i, lT,m, lT,oMt]-triangle tower the [3, 1, l]-simplices focus light 
rays in the angular directions. This is similar to the behaviour of the [1,1,2]- 
and [2, 1, l]-tetrahedra in 2+1 dimensions. Note also that the contributions from 
the terms in square brackets in (j4H) and (j42p are small, because the (absolute 
value of the) prefactor is small (< 0.016 for a > 7/12, which is the range usually 
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considered j3]). This is due to the fact that the extrinsic curvature terms coming 
from the [3, 1]- and the [2, 2]-boundary simplex have opposite signs. 

4. Dynamical triangulation of a black hole 

In this section we will explain how to construct a black-hole geometry in the 
formulation of causal dynamical triangulations. This is not only a good exercise in 
translating metric data (which still depend on a coordinate choice) into invariant 
geometric data, but may also be helpful in deciding on the boundary conditions 
for a path integral over black-hole geometries. It should be kept in mind that 
unlike in Regge calculus |2Z1 "we are working with simplices of fixed squared edge 
lengths and — aa^ for the spacelike and timelike edges respectively. As is 
well known, these fixed building blocks are not well suited for approximating 
arbitrary smooth geometries pointwise, as may be desirable in a numerical study 
of classical Einstein gravity. Instead, as already discussed in the Introduction, 
the application we have in mind here is a nonperturbative quantum superposition 
of spacetime geometries, with classical properties emerging only on sufficiently 
coarse-grained scales. 

This property of dynamical triangulations prevents us from approximating 
a given smooth manifold exactly, even in the limit as the lattice spacing a 
0. For example, it is not possible to obtain a tesselation of three-dimensional 
flat euclidean space using equilateral tetrahedra. Nevertheless, one can usually 
arrange that geometric quantities - for example, the curvature scalar - when 
integrated or averaged over sufficiently large patches match certain prescribed 
values^, and this is what we will employ in the following. 

The explicit construction of a triangulated black hole can be simplified greatly 
by an appropriate choice of coordinates in the continuum. Since causal dynamical 
triangulations come with a distinguished notion of proper time, measuring the 
invariant distance between hypersurfaces of constant time, it is natural to start 
with coordinates in a proper-time gauge, where the lapse function is equal to one 
everywhere. In the following subsections we will describe the relevant black-hole 
configuration (given by the Kottler solution), identify appropriate coordinates 
for it, and construct a triangulation for the r-t-plane of this geometry. The 
latter serves as base space for a full four-dimensional triangulation which can 
be constructed as a triangulation of product type. By its very construction our 
simplicial rendering of the black-hole geometry is highly nonunique, and should 
merely be regarded as an illustration of how it can be achieved. 

^This simple strategy does not work for quantities which are non-negative, for example, 
the square of the Riemann curvature scalar. In such cases, one may need to perform a 
finite "renormalization" or use a more sophisticated way of averaging. Ultimately any such 
prescription must be motivated by physical considerations. 
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4.1 The Kottler solution 



Since tlie metliod of dynamical triangulations requires a positive (bare) cosmolog- 
ical constant, the relevant spherically symmetric black- hole configuration is not 
the Schwarzschild solution, but the so-called Kottler solution which describes a 
black hole in a de-Sitter background with a positive cosmological constant A. The 
line element of the Kottler solution in Schwarzschild-like coordinates {t, r, 6, 0) is 
given by [2H1 

ds' = - [l - — - j^,) dt' + ^1- — - dr' + r'dn' 

=: -f{r)dt^ + firy^dr^ + r^dQ\ (43) 

where M is a mass parameter, L a length parameter defined by = I-, and dQ"^ 
the volume element of the unit two- sphere. Apart from the physical singularity 
at r = there are two coordinate singularities at the positive real roots of /(r) 
(if M < L/a/27, as we will assume from now on), corresponding to the black- 
hole horizon and the cosmological horizon of the de-Sitter background. The 
latter is positioned at r = L if the mass M vanishes. The maximal analytical 
continuation of the spacetime geometry (j43p includes infinitely many black-hole 
and cosmological horizons, but for our present purposes we will only be interested 
in a region inside the cosmological horizon which contains a single black hole. 

To facilitate the translation of the metric data of the Kottler solution to a 
dynamical triangulation, it is convenient to work in a proper-time gauge, that is, 
use coordinates in which the lapse function is equal to 1. A particular set of such 
coordinates are the Painleve-Gullstrand (PG) coordinates (r, r, 6',(^) [2211201, in 
terms of which the metric ()43p assumes the form 

ds^ = -dr^ + (^dr + Vl-/(r)c/r)^ + r^dQ^. (44) 

We observe that the lapse function is equal to 1, but the shift A',, in the radial 
direction is nonvanishing, Nr = a/1 — /. One can also find coordinates with 
the same constant-time surfaces as in (jl^ and A^ = 0, but it turns out that 
for our discussion keeping the nonvanishing shift function is more convenient. A 
remarkable property of the PG-coordinates is the fact that the induced constant- 
time surfaces E,- are fiat Euclidean three-spaces. Our triangulated black hole will 
have the form of constant-time surfaces represented by three-dimensional spatial 
triangulations which are (approximately) fiat and connected to each other using 
the information contained in the shift function Nr. 
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4.2 Triangulation of the base manifold 

Next, we will triangulate the geometry induced on the hypersurface {6 =const., 
(j) =const.}, which is parametrized by a coordinate pair (t, r). The radial co- 
ordinate r is a "proper-length" coordinate in the same sense in which r is a 
proper-time coordinate. It is therefore natural to identify (a rescaled version of) 
r with the discrete time T inherent in a causal dynamical triangulation and (a 
rescaled version of) r with a discrete coordinate R along the one-dimensional 
spacelike triangulation E^. We can introduce a discrete radial coordinate R sim- 
ply by taking a triangulated half-line^ and enumerating its vertices by 0, 1, 2, 
starting from the left-most vertex, say. Since spatial edges by definition have a 
geodesic length a, discrete and continuous radii are related hj r = a R. Similarly, 
one finds r = a^ya + l/AT for the relation between the discrete and continuous 
times. 

Consider now a spacelike edge e in St whose centre has radial coordinate Rt- 
If we want to erect a [2t, lT+i]-triangle over this edge we must decide which vertex 
in the next constant-time line the tip of the triangle should be connected 

to. Because the vector n pointing from the middle of the edge e to the tip of 
the triangle is normal to the hypersurface S^, the displacement of this top vertex 
along Et+1 from the position Rt on is determined by the shift Nr. Explicitly, 
the position Rt+i of the top vertex is given by 

R^^^ = Rt- Nr.{aRT)a-^AT = Rt - Nr{a Rt) ^/a + T/A, (45) 

where At is the modulus of the length of the vector n. In order to construct 
an explicit triangulation, ()45|1 must be approximated by integers. Whichever 
prescription one chooses, the following cases may occur when calculating the 
locations of the tips of [2, l]-triangles: 

(0) The normal vectors ni and n2 of two neighbouring [2, l]-triangles cross 
within the slice AT = 1. For this to happen, drNr must be positive and 
drNr > 1/At i.e. the geometry must vary on scales smaller than a. For the 
Kottler solution this does not occur as long as a -C L = ^/3JA. 

(1) The tips of two neighbouring [2, l]-triangles coincide. In this case the (com- 
ponent in r-direction of the) extrinsic curvature is positive since the normals 
converge toward each other. 

(ii) The tips of two neighbouring [2, l]-triangles with base in are one or 
more edges apart in Sr+i. In this case, we fill the space between the two 
[2t, lT+i]-triangles with the appropriate number of [1^, 2T_|_i]-triangles. The 

^In numerical simulations, there will always be a maximal radius Rmax, because of the 
finiteness of the spacetime volume. 
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extrinsic curvature is zero when a single upside-down triangle is inserted, 
and negative otherwise. 



(iii) The position Rt+i assumes negative values. For the Kottler solution this 
happens whenever the discrete radial coordinate Rt satisfies 

where the approximation in the second step is justified because of a -C L. 
In this case we simply omit the [2, l]-triangles. This way one obtains a 
spacetime boundary which is effectively spacelike (because there are far 
more spacelike than timelike edges), as illustrated by Fig. ^2 It corresponds 
to the spacelike singularity of the Kottler solution at r = 0. 



Figure 11: Example of a dynamically triangulated black hole. The thick lines corre- 
spond to light rays, and the vertical light ray on the right runs along the horizon. 



Fig. ^2 shows the typical features of a dynamically triangulated black hole: ap- 
proaching smaller radii, the tips of the [2, l]-triangles get dragged more and more 
toward the centre of the black hole (left in the figure). Similarly, a freely falling 
particle which starts at some radius r on a surface and normal to it, will 
roughly follow the direction of the timelike normals of these triangles and finally 
fall into the singularity. The horizon in this picture is located at the radius th 
at which the line r = const, is lightlike. 

4.3 Triangulation of the four-dimensional manifold 

Finally, we will sketch how to construct a four- dimensional dynamically triangu- 
lated black-hole geometry by suitably triangulating the fibres F over one of the 
triangulated base spaces B of the previous section. In other words, one has to 
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provide a triangulation of the triangle towers a x F for all triangles a of the base 
space, as described in Sec. 2.3. Since there is no triangulation which is exactly 
spherically symmetric, the spherical symmetry of the Kottler solution can be im- 
plemented at most in an averaged sense. We will formulate conditions on the 
numbers of various simplex types in the fibres which must be satisfied in order 
that the fibre triangulations can be made homogeneous on large scales / ^ a. We 
will not spell out how these building blocks should be distributed over the fibres 
in a maximally uniform manner. 

Since the spatial hypersurfaces S,- of the metric ()44|) are flat, the radii r 
provide a direct measure of the areas of the associated two-spheres {t= const., 
r = const.}. This implies that a vertex tower with radial coordinate R should 
contain (approximately) ciR^ triangles, where ci ~ Att/A^ is of order 1, and 
Aa is the area of a spatial triangle in units of a. This consideration fixes the 
number of triangles in a vertex tower and therefore the numbers of [Sr, Ir+i]- 
and [In, 3R+i]-tetrahedra in a spatial edge tower located between the radii R and 
-R+1. What remains to be specified is the number of [2/j, 2ii;+i]-tetrahedra in the 
spatial edge tower. This can be done by demanding that the (absolute value of 
the) integrated intrinsic curvature of the edge tower be minimal, i.e. as close to 
zero as possible. 

Fixing thus the numbers of simplices appearing in the triangulated surfaces 
of constant time implies that also the numbers of [4^, ly+i]- and [1t,4j'+i]- 
simplices (in the [2t, 1t+i]- and [It, 2T+i]-triangle towers respectively) of the four- 
dimensional triangulation are fixed. Similarly, the number of [Ir.m, ^T,out, 3t+i]- 
simplices in a [2t, lT+i]-triangle tower is already specified by the volume of the 
vertex tower over the tip of the triangle vertex with time coordinate T + 1. All 
we are left with are the numbers of [2T,m, ^T,out, 2t+i]- and [Ir.m, ^r^out, 2t+i]- 
simplices or, equivalently, the numbers of [2^, 2T'+i]-simplices in the timelike edge 
towers [It, 1t+i]- These numbers can also be determined by demanding a van- 
ishing integrated four-dimensional intrinsic curvature scalar, see [12] for details. 

5. Summary and outlook 

In this paper we have derived an explicit expression for the expansion rate of light 
rays for a simplicial manifold of the kind that occurs in the sum over geometries 
in the causal dynamical triangulations approach to quantum gravity. Its prime 
intended use is in a horizon finder in the quantum theory. The vanishing of the 
integrated version of the expansion rate is an indicator for an apparent horizon 
in situations where the geometry along the angular directions is homogeneous 
at sufficiently large scales. Our non-integrated expression for the expansion can 
also be used in more general situations, but a numerical horizon search would 
be considerably more involved, because the expansion rate would have to be 
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monitored not just as a function of radius and time, but also of the two angular 
directions. 

As a first step towards understanding the role of black holes in nonperturbative 
quantum gravity we suggest to investigate the situation where the geometric 
configurations are approximately spherically symmetric. In practice, one would 
first search for the formation of an apparent horizon as function of the boundary 
conditions. This will involve monitoring the integrated expansion rate Sjj^{S) 
given in eqs. pTIjl and as a function of the discrete radius and time coordinates 
of the triangulated base manifold. One will have to work out how the local 
Monte Carlo moves j3] change the counting variables appearing in Sj^{S) in a 
small neighbourhood of triangle towers. In case one manages to find evidence for 
an apparent horizon, one will try to understand whether also other large-scale 
properties of the geometry match those of a (classical) black hole. Ultimately, 
one is of course interested in the quantum deviations from the classical geometry, 
especially near the horizon and for very small radii, to obtain further insights into 
a possible quantum origin of (horizon) entropy and a quantum resolution of the 
central singularity. These are doubtless ambitious goals, but with some luck and 
ingenuity they may just be within reach of our computational means. 
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Appendix 

In this appendix, we define the concept of affine coordinates for flat simplices, and 
use them to describe hypersurfaces of constant time, as well as the surfaces S 
used in the computation of the extrinsic curvatures in Sec. 3.3. For illustrative 
purposes, we also perform a simple model calculation for d = 2. 

For the computation of various geometrical quantities like the angles ()34|1 
and appearing in the expansion rate (jHSl), it is very convenient to introduce 
affine or barycentric coordinates (see, for instance, |SI])- Consider a (i-simplex a'^ 
embedded into (i- dimensional Minkowski space. Let the vectors Vj, j = 1, . . . , d+1 
point from the (arbitrarily chosen) origin to the (d+l) vertices of the simplex. 
One can then describe an arbitrary point P in cr'^ as the centre of mass of (d+l) 
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appropriately chosen masses > located at the vertices of the simplex, 

d+l 

P = Y^m^Vj, (47) 
i=i 

where the sum of the coordinates is normalized to 

^m^ = l. (48) 
i=i 

An affine vector can be defined by the difference of two points Pi and P2, 

d+l 

= J2 {mi - m{)vj . (49) 
i=i 

The sum of the vector components := {ml — m{) in ()49p vanishes because of 
the normalization ()48|) . We can therefore replace the vectors vj in fl49p by the 
affine basis vectors 

d+l 



k=l 

which are overcomplete due to Ylj^j — ^- 

In order to express the metric in affine coordinates we still need a basis of 
one-forms dual to {aj}. Because of the overcompleteness of the vector basis, the 
dual basis {a-'} is defined by 

(a,,a'=)=5j:=5j-^, (51) 

where (5* is the projector onto the space of affine coordinates, i.e. it projects an 
arbitrary ((i+l)-tuple of numbers into one whose sum vanishes and it acts as the 
identity on tuples (?/^, . . . , y'^^^) which have a vanishing sum ^ y^ = 0. 

The (Minkowski) metric components of a simplex in affine coordinates take 
the form [HI] 



1 



where /f • is the squared geodesic length of the edge between the i-th and j-th 



vertex. The metric two- form is then given by r/ = T^y-a^a-'. Since the (5*'s on 
the right-hand side of are projections onto the space of affine coordinates 
they can be replaced by ordinary Kronecker symbols 5j if fl52p is contracted with 
coordinate tuples that already fulfill the affine coordinate condition (i.e. sum 
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up to zero). One can check eq. ()52|) by contracting it with the simplex edges 
e(jj) := Vj — Vi = SLj — aj which gives the expected resuh 

Vie(ij),e^ij)) = 1^. (53) 

With the help of the barycentric coordinates one can give an easy characteri- 
zation of the fibres with respect to the two different product structures introduced 
in Sec. 2, namely, the constant-time hypersurfaces and the surfaces S of codi- 
mension two. First, consider a (i-simplex of type [Nt, Nt+i], that is, a simplex 
having Nt vertices in the hypersurface and Nt+i vertices in the neighbouring 
surface T^t+i- The barycentric coordinates (m^, . . . , m'^"'"^) of any point in the 
intersection of T^T+t, < t < 1, and the d-simplex [Nt, A^r+i] satisfy 

Nt d+1 

^ m^' = 1 - t = 1 - m?. (54) 

i=l j=ArT+l 

This condition describes a linear subspace whose points have a "constant dis- 
tance" from the T- and the (T+1)- vertices. 

A surface 5* can be characterized in an analogous manner. Consider a d- 
simplex of type \NT,in,^T,out,^T+^\, with A^T,m vertices in the inner and Nt^oui 
in the outer vertex tower at time T and Nt+\ vertices in the vertex tower at time 
T+1. (We use "inner" and "outer" merely as labels to distinguish between the 
two vertex towers at equal time.) Then the intersection of the d-simplex with 
the surface S is the set of points having "constant distance" to these three set of 
vertices, 

^T,i-n NT,in+NT,out d+l 

= r, ^ = s, = t, (55) 

where r + s + t = 1. Obviously these points also fulfil eq. (j5l|l . which shows that 
S* is a sub manifold of T^T+t- 

To illustrate the construction, consider the simple case of a two-dimensional 
Lorentzian triangulation containing a spacelike hypersurface T^T+t- We want to 
calculate the angles p defined in (jSH) to compute the extrinsic curvature of Sr+t- 
Take a [2^^, ly+ij-triangle and introduce barycentric coordinates {im} ,m? ,m?). 
From eq. (|52|) . the components of the Minkowski metric in affine coordinates are 
given by 

{0 for z = j, 
— Y for e(jj) a spacelike edge, (56) 
^ for e{ij) a timelike edge, 

where we have assumed that the metric is contracted with affine coordinate tuples 
only. The tangential vector to T^T+t is s = c<j(l,— 1,0), with a normalization 



39 



constant Cg. Orthogonal to this and future-directed is the normal vector n = 
Cn(— |, — |, !)• If we take as the timelike boundary simplex the edge 6(1^3), the 
future-directed vector parallel to it is eg = Co(— 1,0, 1). Lastly, orthogonal to Bq 
and pointing into the 2-simplex is the vector 

e. = c.(^— ,1,-). (57) 

Application of the formulas ()34j) to the (normalized) vectors Bq, §1 and n leads 
to the angle 

P[2,i] = 

By an analogous calculation for an "upside-down" triangle [l-r, 2^-1-1] one obtains 
the corresponding angle p[i^2] = ~P[2,i]- 
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